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Abstract 

In the present paper we study the oscillation spectrum of neutron stars containing both 
ordinary matter and dark energy in different proportions. Within the model we consider, 
the equilibrium configurations are numerically constructed and the results show that the 
properties of the mixed neuron-dark-energy star can differ significantly when the amount of 
dark energy in the stars is varied. The oscillations of the mixed neuron- dark-energy stars 
are studied in the Cowling approximation. As a result we find that the frequencies of the 
fundamental mode and the higher overtones are strongly affected by the dark energy con- 
tent. This can be used in the future to detect the presence of dark energy in the neutron 
stars and to constrain the dark-energy models. 
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1 Introduction 



The cosmological observations of the present Universe provide evidences for the existence 
of a mysterious kind of matter, called dark energy, which governs the accelerated expansion 
of the Universe [1], [2]. The dark energy constitutes 73% of total energy content of the 
Universe and it exhibits some unusual properties such as negative pressure to density ratio 
w (in hydrodynamical language). 

The fundamental role that the dark energy plays in cosmology naturally makes us search 
for local astrophysical manifestation of it. If dark energy was discovered at astrophysical 
scales it would give us new tools for more profound investigation of it and its properties. 
For example, the existence of dark energy makes us expect that some of the present-day 
stars are a mixture of both ordinary matter and dark energy in different proportions. The 
study of such mixed objects is a new interesting problem in the current investigations, see 
for example [3]- [11] and references therein. 

The natural question that arises is how we can detect the presence of dark energy in 
the stars. Are there any imprints of dark energy that can be observed and how? In present 
paper we consider one such possibility based on the gravitational waves radiation by the 
neutron stars. It is well known that the characteristics, and more precisely the spectrum of 
the gravitational waves emitted by the compact objects, depend on the interior structure 
of the compact objects. This is the base of the asteroseismology - the interior structure 
of the compact objects can be revealed by investigating the spectrum of the gravitational 
waves they emit [121 [13]. Currently a lot of efforts are devoted to the detection of the 
gravitational waves. Several ground-based detectors as LIGO, VIGRO, TAMA300 and 
GEO600 are operating and the launching of space-based detectors, for example LISA, are 
expected in near future. The detection of the gravitational waves will open a new window 
to the Universe and will provide us with a unprecedented tool for studying the compact 
astrophysical objects. 

In the present paper we consider a certain model of mixed neutron-dark-energy stars 
(MNDES) and study their oscillation spectrum in the Cowling approximation. We show 
that the spectrum of MNDES contains imprints of dark energy and could be used to detect 
the presence of dark energy in the neutron stars. 

2 Equilibrium model for mixed neutron-dark-energy 
stars 

The nature of the dark energy is a mystery at present and it is a great challenge to current 
physics to solve this problem. Nevertheless, in their attempts to understand and model 
the dark energy, theoretical physicists have adopted some effective descriptions of it. The 
most simple and direct description is by a perfect fluid with a negative pressure. The other 
effective descriptions of the dark energy depend on the value of w. If w > — 1 a possible 
description is provided by a scalar field with an appropriately chosen potential. In the 
case when w = —1 the dark energy can be explained by the cosmological constant. If 
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w < —1 then a possible description of dark energy is provided by scalar fields with negative 
kinetic energy, the so-called phantom (ghost) scalars [13], [U]. The negative kinetic energy, 
however, leads to severe quantum instabilities and this is a formidable challenge to the 
theory. However, there are claims that these instabilities can be avoided [16]. In general, 
the problem could be avoided if we consider the phantom scalars as an effective field theory 
resulting from some kind of fundamental theory with a positive energy [T7], [IB]. In this case 
the phantom scalar description of the dark energy is physically acceptable. In this context it 
is worth noting that the phantom-type fields arise in string theories and supergravity [21] 

At present it is not clear whether w < —1, w = —1 or w > —1. It seems, however, that 
the current observational constraints favor w < —1 j25]^ [30]- That is why in the present 
work we adopt an effective description of the dark energy by a phantom scalar. Since a very 
little is known about the interaction of dark energy with the normal matter we consider a 
simple model with minimal interactions for the phantom field - we do not include terms 
describing non-minimal interaction between the phantom field and the normal matter in 
the field equations [11]. Under these conditions the Einstein equations in the presence of 
dark energy read [TT] 

R^i, = 8n{T^^ - ^Tg^^) - 2d^ipd^ip, (1) 

Here T^i, is the energy-momentum tensor of the ordinary matter in the perfect fiuid descrip- 
tion with energy density p and pressure p: 

T^u = {p + p)u^u^ + pg^^. (2) 

For the ordinary matter we impose the natural conditions p > p > 0. The density of the 
dark energy sources is denoted hj po- The dark energy sources will be called dark charges. 
When we consider the local manifestation of dark energy on astrophysical scales, i.e. scales 
much smaller than the cosmo logical scales, the phantom potential V({(p) can be neglected 
and that is why we have set U{(p) = 0. 

The contracted Bianchi identity applied to the field equations gives 

v.t; = pdV^^. (3) 

This equation projected orthogonally to the 4- velocity u'^ gives the equation of motion of the 
fiuid. In the case of static configurations (not necessary spherically symmetric) it reduces 
to the equation describing the hydrostatic equihbrium: 

{p + p)di\n^/\g^\+ dip = podiif, (4) 

^The perfect fluid description of tlie dark energy also suffers from instabilities due to the imaginary 
velocity of the sound. From a classical point of view the massless phantom field is even more stable than 
its usual counterpart [19], [20] . 
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where the subscript i is for the spatial coordinates. For barotropic equation of state (EOS) 
of the ordinary matter p = p{p), eq. (jH]) can be rewritten in the form 



d\n^\+d [ ^ = -P^dip. (5) 
J p + p p + p 

The integrabihty condition for this equation gives 

PD = /M(p + P), (6) 

where f{^p) is arbitrary function of ip. In the present paper we consider the simplest choice 
for /(v?), namely f{<p) = k where > is a constant, i. e. 

pn = k{p + p). (7) 

In fact the constant k is a measure for the content of dark energy in the star, while 
(p + p) determines the dark energy profile in the star. 

Now let us consider the projection of ([3]) along the 4-velocity u^. We find 

-|-(p + p)VX = Pz,^. (8) 

where ^ is the derivative along the fluid 4-velocity. Combining this equation with the 
baryon conservation law {nu'^) = we find 



_dp^ (p + p) dn ^ rf^ 
dr n dr dr^ 



or equivalently 



dp = ~^ dn — podip, (10) 
n 

which is the thermodynamics first law for our system of ordinary matter and dark energy. 
Hence we conclude that 



p = p{n,(p), (11) 

i.e. the energy density is a function not only of the particle number density but depends 
also on the scalar field. In order to find the function p{n, (p) we use the barotropic equation 
p = p{p) and eq. (171) which substituted in ( JTOl) and some algebra gives 
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dp 
P + P 



dn^ 



(12) 



where 



(13) 



Eq. (fT2!) determines the function p{n, ip) and it is easy to see that if in the absence of dark 
energy p{n) = F{n), then p{n, (f) = F{n^) = F{e~'"^n) in the presence of dark energy. 

From now on we will focus on spherically symmetric equilibrium solutions. In this case 
the spacetime metric can be written in the well-known form 



ds^ = -e^^(^Ut^ + e^^^'-^rfr^ + r\de^ + sin' Od^'') 
and the reduced field equations are the following 



(14) 



[r(l-e--)]=8.pr^-r^ (^)' 



dr 

d$ _ 1 
dr 2r 

d 
dr 



-2A 



e'^ -1 + Snpr^e 



2^2A ^2 I 



dr I 



(15) 



dr 



Anpor e 
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together with the equation describing hydrostatic equilibrium @ in the spherically sym- 
metric case 



, (i$ dp du) 
dr dr dr 



(16) 



Here it is worth noting that the system of equations f|T5l) - f|T6l) . can be viewed as de- 
scribing an equilibrium spherically symmetric configuration with an effective anisotropic 
pressur^. As a hole we can say that the presence of dark energy in the star yields a 
pressure anisotropy of certain kind and can be considered as an effective quasi-local 
equation of state for the anisotropic part of the pressure. The local anisotropy also affects 
the thermodynamics first law as we have already seen in our case. 

We numerically solve the system (fT5l) - (fT6l) with the following boundary conditions 



A(0) = 0, p(0) = po, ^(0) = 0, ^{oo) = 0, $(oo) = 0. 



(17) 



For a discussion of locally anisotropic selfgravitating systems we refer the reader to [21] ■ 
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The radius R of the star is determined form the condition p{R) = 0. We also use a 
polytropic equations of state for the ordinary matter 



P = p{n^) = Knouib ( — ] = Knonib i e ] , (18) 

p = p{n^) = n^nib + _ ^ = e'^'^nnii, + ^ _ ^ , (19) 

where n is the particle number density, = 1.66 x 10~^*^(y', no = 0.1/m~'^. We have chosen 
r = 2.46 and K = 0.00936 obtained when fitting the tabulated data for EOS II [32j, but 
the results are qualitatively the same for other values K and F. 

The mass, the number of baryons and the dark charge of the star are given by [11] 



M 



--T I RlV^d^x = [ (P + ^p)y/^d^x = A7i [ (p + 3p) e^+^r^rfr, (20) 

471" J star J Star Jo 



R 



N = Att ne\'dr, (21) 
Jo 

D = [ pDV^d^x = An [ pDe^+\'^dr = — I V^ifdT.^'. (22) 

J star Jo 47r J 3-2^ 

Other equivalent formulas for calculating M and D can be obtained by using eqs. f llSI) . 

The numerical results for the background solutions are shown on Fig. [1] where the mass 
M as a function of the central density po and of the radius R are presented. On Fig. |2]the 
dark charge D and the ratio of the dark charge to the total mass D/M of the MNDES as 
functions of the central density are shown. The sequences of solutions are obtained when 
varying the value of the cental density of the ordinary matter po while keeping the dark 
energy parameter k fixed. As we can see the properties of the star can vary significantly 
when the portion of the dark energy inside the star increases, i.e. when we increase k. 
For large values of k (i.e. large portions of dark energy) the total mass of the star can 
be much larger than the pure neutron star mass for the same equation of state and the 
radius of the star also increases. These effects are due to the phantom field which yields 
repulsive rather than attractive force. In fact the additional force on the right hand side of 
eq. (fT6|l is repulsive, Pd^ > and this can be easily seen form the field equation for ip. In 
hydrodynamical language, the additional force due to the anisotropic nature of the fluid is 
pointing outward. We have found numerically that there are no equilibrium configurations 
for k >1 since the repulsive force dominates on the attractive gravity in that case. 

Now let us turn to the question about the stability of the equilibrium configurations. The 
dynamical stability is usually studied by considering their radial perturbations and reducing 
the problem to an eigenvalue problem. Since the equilibrium configurations in our case are 
parameterized by one parameter (for fixed k), say the central density of the ordinary matter 
Po, the eigenvalues will be continuous functions of them. What is important in the present 
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Figure 1: The mass of the MNDES as a function of the central density (left panel) and of 
radius of the ordinary matter component (right panel) for several values of the dark energy 
parameter k {k = corresponds to a pure neutron star). 



context is the point between the stable and instable configurations. On that point we must 
have a zero eigenvalue or in other words a static perturbation transforming the equilibrium 
state with po to the equilibrium state with po + ^Po- The static perturbation preserves the 
number of particles N of the ordinary matter. Therefore the perturbed equilibrium state 
have the same N as the initial equilibrium state which means that 



dN 

6N = —6po = 0. (23) 
dpo 

Therefore the change of stability is where |^ = 0. Moreover, for our model and for fixed 
k, one can show thati 

6M = j26N, (24) 

where fi is the red-shifted chemical potential, p = ^^e*. It can be shown from ( !T6|) that 
p, is constant throughout the star. Consequently the mass and the particle number peak 
at the same location, as one can see form Fig. [3l Therefore the stability change occurs at 
the critical energy density pg for which |^ = 0. The configurations with po < Pg are stable 
while those with po > Pg are instable because the configurations with pg < pg have lower 
energy than those with pg > pg as one sees on Fig. [31 



3 Perturbation equations in Cowling approximation 

In this section we derive the equations describing the non-radial perturbations of the MN- 
DES in the so-called Cowling approximation [33]- [31]. In Cowling approximation the 

''More precisely we have 5AI — fiSN + DSipao-, however we consider fixed cosmological value (poo- 
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Figure 2: The dark charge in solar mass units (left panel) and the ratio of the dark charge to 
the total mass of the MNDES (right panel) as functions of the central density. The results 
are for the same solutions as those shown on Fig. [H 




Figure 3: The dependence of the star mass on the ordinary matter particle number for 
different values of k. 
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non-fluid degrees of freedom are kept flxed - for example, within the framework of general 
relativity, the spacetime metric is kept flxed. Despite of this simpliflcation the Cowling 
formalism turns out to be accurate enough and reproduces the oscillation spectrum with 
good accuracy. In fact the comparison of the oscillation frequencies obtained by a fully 
general relativistic numerical approach and by the Cowling approximation shows that the 
discrepancy is less than 20% for the typical stellar models [55]. The Cowling approxima- 
tion is trustable even in the case of presence of scalar flelds [SB]. We apply the Cowling 
approximation to our model by keeping the metric and the scalar phantom fleld frozen. 

The equations describing the perturbations in Cowling formalism are obtained by varying 
the equation for the conservation of the energy-momentum tensor ([3]). Taking into account 
that the metric and scalar fleld are frozen, we flnd V^5T^ = where 

<5T; = {5p + 5p) u^v" + (p + p) {u^bu' + bu^u') . (25) 
Projecting the equation VuSTj^ = along the background 4-velocity we have 

u'VJp + V, [(p + p)6u''] + {p + p) aju-" = 0. (26) 

Projecting orthogonally to the background 4-velocity by using the operator -P^ = 5^ + 
u'^Ufj,, we obtain 

{6p + 6p) a^ + {p + p) u" {VM^ - V^5u^) + V^6p + u^u^'V^p = 0, (27) 

where = u'^'VuU^ is the background 4-acceleration. Taking into account that u = 
(u*, 0, 0, 0) the above equation can be rewritten in the form 

(Sp + Sp) ai + {p + p) u^dtSui + di5p = (28) 

with z = 1, 2, 3 = r, 6*, 

At this stage we can express the perturbations of the 4-velocity via the Lagrangian 
displacement vector namely 

dt M* 

Now let us consider eq.( 125]) for i = 1 and i = 2. Since ag = = we flnd 

{p + p){u')^dUe + de6p = 0, (30) 

(p + p)(u'yd^U + d^Sp = 0. (31) 
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(29) 



Taking into account that p,p and -u* depend on r only, the integrabihty condition for the 
above equations gives 



(32) 



From this condition and the fact that the background is spherically symmetric we find that 
^0 and ^(f, are of the form 



Im 

^4> = -J2^iUr,mYi 



(33) 
(34) 



Im 



where Yim{6, (p) are the spherical harmonics. From now on, in order to simplify the notations, 
we will just write ^ = —VYim when we have expansion in spherical harmonics. 

We proceed further with finding the expressions for the density and pressure perturba- 
tions. From eq.f l27p we have 



u^5p 



-9 



(35) 



It is convenient to express in the form 



r = e-^—Yim 



,2 



(36) 



and substituting in the above equations, after some algebra we fincS 



Sp= -(p + p) 



dr 



dr 



Imi 



(37) 



where in the last step we have taken into account that = and eq. f|T6l) . 

In order to find the perturbation of the pressure we first use the relation between the 
Eulerian and Lagrangian variations, namely 



Sp = Ap- CdrP (38) 
with being the Lagrangian variation. From the equation of state we have 



^The derivative with respect to the radial coordinate r wiU be denoted by prime or by the standard 
symbol interchangeably. 
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(39) 



In this way we obtain the following formula for the perturbation of pressure 



.W\ dp .W 

ar ) ar 



(40) 



Having the exphcit expressions for 6p and 6p we can put them into eg. (1281) and after 
lengthy calculations we obtain the dynamical equations for W and V 



,A-2<I> 



{p + p)—^d^W 



d ( , .dp 
ar { dp 



dp _pW dp dip _p^W 
dr dp^^ dr 



ar ar 



W 



\ dp J dr 



(41) 



{p + p)e~'''dtV = (p + p) 



dp 
dp 



ar 



W dp d<p _^ 
r-^ dp dr r 



W 

~~2' 



(42) 



From now on we will assume for the perturbation functions a harmonic dependence on 
time, i.e. W{r,t) = l^(r)e''^* and V{r,t) = V{r)e^'^K The system (I3I])-(112]) can be consider- 
ably simplified by combining the equations in an appropriate manner. When differentiating 
equation fH2|) and adding it to equation fj4T]) . and also using eq.f lT6|) . we find 



V = 2$V - e^^ - '-^pn^V. 

r'' P + P dr 



(43) 



This equation together with the equation (1421) solved for W, form a system which is equiv- 
alent to (I4ip - (H2]) but much more tractable: 



dp ^ -' p + p dr 



(44) 



dp 



W 1 + 1 d^^^ 

— ^ — PD-r^- 

p + p dr 
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The boundary condition at the star surface is that the Lagrangian perturbation of the 
pressure vanishes 



(P + p) L'^e-^'-V + t'e-"'-^ - -^%-''%) = 0. (45) 

y T'^ p + p dr J 

The boundary conditions at the star center can be obtained by examining the behaviour 
in the vicinity of r = 0. For this purpose it is convenient to introduce the new functions W 
and V defined by 

W = Wr^+\ V = VrK (46) 
Then one can show that at r = the following boundary condition is satisfied 

W = -IV. (47) 

4 Numerical results for the oscillation spectrum 

As we have shown in the previous section the properties of the MNDES depends strongly 
on the dark energy content of the star and the masses and the radii can differ significantly. 
It is natural to expect that this will influence the oscillations of the star. In order to find 
the oscillation spectrum of MNDES we have numerically solved the spectral problem given 
by the differential equations (jSj) together with the boundary conditions (H5l) and ( H71) . In 
this section we will present the results for the /-modes of the MNDES and for the higher 
fluid modes - pi and p2- AH the dependencies are shown up to the maximum mass for the 
corresponding values of the parameters, i.e. where a change of stability is observed. 

For pure neutron star an empirical relation between the average density and the /-mode 
frequency exist^ which does not depend much on the EOS |T2]. As we can see from Fig. H] 
this relation differ for MNDES with different dark energy content (i.e. for different values 
of fc), but it remains relatively close to the pure neutron star (the case with = 0). Only 
for large values of k the relation changes more significantly - the frequencies grow slower as 
a function of the average density. 

The influence of the dark energy on the oscillation spectrum is more pronounced on the 
Figs. El and [7] where it is shown the oscillation frequencies / and the normalized frequency 
cu as a function of the total mass M of the MNDES for the /, pi and p2 modes. The values 
of the frequencies / decrease when we increase the dark energy content for all of the fluid 
modes and the differences with the pure neutron star case {k = 0) can be large. This effect 
is more pronounced for the higher fluid modes and for larger values of k the frequencies of 
the p2-niode can reach two times the frequencies of the pure neutron star case as we can see 

"The /-mode frequency is proportional to the average density. 
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on Fig. [71 This is a very strong effect and by observing more than one ffuid mode of the 
same neutron star we can put serious constraints on the dark energy content in the MNDES 
and consequently on the dark energy model itself. 

An interesting observation is that the normalized frequency uj behaves differently as a 
function of the total mass M for the different fluid modes as we can see from the right 
panel on Figs. I3 E] and [71 The results show that the normalized u increases when we 
increase the parameter k for the /-mode whereas for the higher fluid modes u decreases 
with the increase of the dark energy content. This effect is qualitatively different from most 
of the other alternative models of neutron stars where the qualitative behaviour of the 
normalized frequencies u is the same for all of the modes when we vary the parameters 
of the corresponding model. This property can be potentially used to distinguish between 
MNDES and other alternative neutron star model. Again the differences between oscillation 
frequencies of the pure neutron star and of the MNDES are larger for the higher order modes. 

These properties of the oscillation spectrum can be easily seen also on Fig. [HI where the 
normalized frequencies of the /, pi and p2-niodes are shown as a function of k for fixed 
value of the MNDES mass M = IAMq. The results show that when we increase the value 
of k not only the absolute values of the frequencies change in a bigger range for the higher 
order fluid modes, but also the relative change in the frequencies is larger. On Fig. [HI we can 
also observe one of the properties of the oscillation spectrum that we discussed above - the 
normalized frequency increases when we increase the parameter k for the /-mode whereas 
for the higher fluid modes u decreases with the increase of k. 

5 Conclusions 

In the present paper a model of neutron star containing dark energy and its oscillation 
spectrum are studied. We model the dark energy part of the star by using a phantom 
scalar field [H], [TB], [TT]. The neutron-dark-energy star solutions are obtained by solving 
numerically the reduced system of field equations. The results show that the properties of 
the MNDES can be significantly different from the pure neutron star case. One of the main 
difference is that the masses and the radii can be significantly larger in the presence of dark 
energy. The stability of the MNDES is also studied and it turns out that similar to the 
pure neutron star case a change of stability is observed at the critical energy density Pq for 
which 1^ = 0. 

The second part of the paper is concentrated on examining the possible dark energy 
imprints in gravitational wave spectrum of the obtained MNDES. In order to do that we 
study the oscillation spectrum of the obtained solutions in the Cowling approximation where 
the non-fluid degrees of freedom are kept fixed. The values of the oscillation frequencies of 
the MNDES can differ significantly from the pure neutron star case and the differences are 
bigger for the higher fluid modes. Thus the observation of gravitational waves emitted from 

**Such alternative models are for example the neutron stars in the scalar-tensor theories [36j or in the 
tensor- vector-scalar theory of gravity [37] . 
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Figure 4: The /-mode frequency as a function of the average density of the ordinary matter 
component \/ MjFt?. The results for several values of the dark energy parameter k are 
shown. 
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Figure 5: The frequency / as a function of the mass M (left panel) and the normalized 
frequency a; as a function of M (right panel) for the / mode. The results for several values 
of the dark energy parameter k are shown. 
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Figure 6: The results for the pi-mode of the same solutions as shown on Fig. [5l 




Figure 7: The results for the p2-niode of the same solutions as shown on Fig. [5l 
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Figure 8: The normalized frequency w as a function of k for fixed value of the mass M = 
IAMq. The results for the /, pi and p2 modes are shown. 

neutron stars could be used to detect the presence of dark energy in the neutron stars and 
to put constraints on the dark energy properties in general. 

We will finish by mentioning a possible further extension of the present work. Since the 
scalar phantom field violates the weak/null energy condition it can maintain configurations 
with nontrivial topology - the so-called wormholes. Such a solution describing a star har- 
bouring a wormhole at its core was recently found in ^U\. It would be very interesting to 
study the influence of the possible nontrivial topological structures in neutron stars, like 
wormholes, on the gravitational wave spectrum of the stars. 
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